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Abstract 

We present comprehensive numerical results for domain growth in 
the two-dimensional Random Bond Ising Model (RBIM) with non- 
conserved Glauber kinetics. We characterize the evolution via the 
domain growth law, and two-time quantities like the autocorrelation 
function and autoresponse function. Our results show evidence for 
the crossover from a pre-asymptotic regime with "power-law growth 
with a disorder-dependent exponent" to an asymptotic regime with 
"logarithmic growth". We compare this behavior with previous re- 
sults on one-dimensional disordered systems and we propose a uni- 
fying picture in a renormalization group framework. We also study 
the corresponding crossover in the scaling functions for the two-time 
quantities. Super-universality is found not to hold. Clear evidence 
supporting the dimensionality dependence of the scaling exponent of 
the autoresponse function is obtained. 
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1 Introduction 



The coarsening dynamics or domain growth of a thermodynamically unstable 
system, such as a ferromagnet or a binary mixture quenched below the critical 
point, is a well-established paradigm of slow relaxation. From the pioneering 
work of Marro et al. [TJ up to the reviews of Bray j2] and Puri [3] , the equal- 
time properties have been extensively investigated. Recently, the interest has 
shifted to the study of two-time quantities, since the ordering process offers a 
simple and non-trivial arena for the understanding of aging phenomena [HE]. 
There is also a long-standing interest in studying the robustness of the basic 
features of coarsening systems when quenched disorder (without frustration) 
is introduced |H UJ- Taking the Ising model with ferromagnetic couplings 
as the basic system, the class of disordered systems investigated includes 
the random field Ising model (RFIM) [SJ and the random bond Ising model 
(RBIM) [7J. 

The first effect of disorder is to change the relaxation mechanism from 
purely energy lowering to activated, due to the creation of energy barri- 
ers [9]. The question is to what extent this affects the growth law, namely 
the dependence on time of the linear size of domains L(t). Intuitively, the 
presence of energy barriers is expected to slow down the coarsening process. 
However, the problem turns out not to be so straightforward, since the bar- 
rier size may depend on the domain size, thus determining the nature of 
the asymptotic growth law through a non-trivial feedback process. So far, 
this issue has been well understood in the case of one- dimensional systems. 
The d = 1 RFIM [TO] and the d = 1 RBIM [llj offer instructive and clear 
instances of disorder forming either L-dependent barriers (in the RFIM) or 
L- independent barriers (in the RBIM). In both cases, there is a crossover. 
In the d — 1 RFIM, it occurs from pure-like power-law growth to disorder- 
dominated logarithmic growth. In the d — 1 RBIM, the crossover takes place 
in the reverse order: from a pre-asymptotic disorder-dependent power law to 
an asymptotic pure-like power-law growth. 

In higher dimensions, a crossover from algebraic to logarithmic behavior 
has been found in the growth law of the RFIM [T2J H3J [HJ [T5] . A crossover 
from algebraic to slower than power-law behavior, possibly logarithmic, has 
been reported also in the case of the random site Ising model (RSIM) [161 EJ 
IT?] . According to the phenomenological theory of Huse and Henley (HH) [18J, 
this means that in these systems energy barriers ought to scale like a power 
of L. 
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Instead, in the d = 2 RBIM the role of barriers and the nature of the 
growth law has remained controversial, since the available numerical evi- 
dence shows a steady algebraic growth with a disorder- dependent dynamical 
exponent [HI [2TJ [22]. The open issue is whether the data manage to probe 
the genuine late time behavior, in which case energy barriers ought to scale 
logarithmically with L |19] , or the HH paradigma is valid also for this sys- 
tem, in which case the data would be probing a long pre-asymptotic regime, 
eventually to be followed by a crossover to logarithmic growth [23| |2"%]. That 
an HH-type scenario might yield an effective algebraic law at intermediate 
times had already been noticed by Bouchaud et al. [25] in the context of spin 
glasses. Our contribution in this paper is precisely the clearing up of this 
point, producing the numerical evidence for the existence of the crossover to 
logarithmic growth in the d = 2 RBIM. 

A related important question is how disorder affects the scaling functions, 
namely whether it enters only through the modified growth law or if there is 
also an explicit disorder dependence. In the first case, scaling functions ought 
to be the same for pure and disordered systems, once time is reparametrized 
through L(t). This is referred to as the super-universality (SU) hypothe- 
sis [26] , whose validity has been supported mainly from numerical results for 
the equal-time correlation function or structure factor [6j [HI [151 IE]- More 
recently, the validity of SU has been extended by Sicilia et al. [20] to the 
geometrical properties of domain structures and by Henkel and Pleimling 
(HP) [2"Tj |2"2"] to two-time quantities, such as the order parameter autocorre- 
lation and autoresponse functions in the d = 2 RBIM. However, next to the 
features supporting SU validity, there exists also some evidence showing SU 
violation, as reported by HP [22] in the short distance behavior of the d = 2 
RBIM and by Park and Pleimling [17J in the space-time correlation function 
in the RSIM. What is missing, here, is the general framework where to fit 
coherently these conflicting pieces of information, with the end-result that 
the issue of SU validity is regarded as yet to be understood. 

By contrast, a clear picture emerges from the d — 1 cases mentioned 
above [TUJ |TT], where the crossover can be investigated with considerable 
precision and turns out not to be limited to the growth law, but to take 
place also in the correlation and response functions. Then, the issue of SU 
validity or lack of validity is reduced to deciding whether the asymptotics are 
or are not disorder-dominated. Adopting this point of view, our finding of 
the crossover in the growth law in d = 2 RBIM strongly suggests that in this 
system too there might be full violation of SU and that this might fit into the 
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pattern extracted from the d = 1 systems. Indeed, this is what we find from 
the autocorrelation and autoresponse function data and, putting together 
these new data with those for d — 1 systems, we have succeded in con- 
structing a new and coherent picture of general character for the coarsening 
dynamics of weakly disordered systems. This is most effectively formulated 
into a renormalization group framework in terms of the competition between 
a pure fixed point and a disordered fixed point. Then, the variety of different 
scenarios observed arises from the different possible choices in the relative 
stability of these fixed points. 

Finally, the presence of disorder allows for an exhaustive investigation of 
the scaling properties of the dynamical susceptibility. The data in the pre- 
asymptotic regime of the d = 2 RBIM, where the growth law is algebraic 
with good precision, together with those for the d = 1 RBIM, produce con- 
clusive evidence for the ci-dependence of the exponent involved. Our results 
are in agreement with the conjecture put forward in studies of the pure sys- 
tem [51 [27J 12H1 129] . A preliminary account of this result has been presented 
in Ref. [II]. 

This paper is organized as follows. In Sec. [2J we provide a general frame- 
work, based on renormalization group ideas, to understand the effects of 
disorder on domain growth. In Sec. EU we summarize our understanding of 
coarsening in d = 1 disordered systems. The bulk of the paper is contained in 
Sec. HI where we present comprehensive numerical results for domain growth 
in the d = 2 RBIM. We conclude the paper with a summary and discussion 
in Sec. [5j 

2 General Framework 

This section is devoted to an overview, using scaling arguments and renor- 
malization group (RG) terminology, of the possible deviations from pure-like 
behavior due to disorder. 

2.1 Growth Law 

The linear size of domains in pure systems asymptotically grows according 
to the power law 

L(t) = J D t 1/2 , (1) 
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which holds when time t is large enough for L(t) to dominate all other lengths 
in the problem. Here, D is a prefactor which is dependent on temperature, 
and z is the dynamical exponent. It is straightforward to verify that the 
above power law satisfies the homogeneity relation 

L(t) = IL (1) (2) 

where I is arbitrary and can be regarded as the scale parameter in an RG 
transformation. In this paper we shall consider only non-conserved dynamics. 
In pure systems the dynamical exponent is z = 2, independent of the system 
dimensionality [2J. 

Assuming scaling, the growth law for disordered systems is obtained gen- 
eralizing Eq. (j2J) to 

= ( 3 ) 

where e is a parameter quantifying the amount of disorder and is the 
corresponding scaling exponent. Once the growth law is written in this form, 
the standard crossover scenario follows from the competition between the 
pure fixed point at e = and the disordered fixed point at e = oo. The 
relative stability of these two fixed points is regulated by the sign of <ft, with 
the pure fixed point being stable or unstable for > or < 0, respectively. 
In the first case the asymptotic growth law would remain the pure one, while, 
in the second case, there would be the crossover to an asymptotic behavior 
different from that of Eq. fl2J and controlled by the disordered fixed point. 

In place of investigating directly the sign of 0, it turns out to be more 
convenient to introduce the crossover length A(e) ~ e 1 ^ and to rewrite Eq. 

l(M)=K'(^f). (4) 

At fixed points the crossover length values are either A* = or A* = oo, and 
since A(e) flows from oo toward 0, the fixed point corresponding to A* = is 
stable, while the one with A* = oo is unstable. Therefore, the relevancy or 
irrelevancy of disorder can be deduced from the behavior of A(e) as e — > 0. 
Specifically, since e = corresponds to the pure fixed point and the two 

1 As we shall see, in the case of the d = 1 RBIM Eq. ((4]) applies also to the more general 
case where A(e) does not scale like a power of e. 
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possible limits of A(e) are 



limA(e) = \ (5) 

if the first one applies, the pure fixed point is stable, while if the second limit 
applies, the disordered fixed point is stable. 

As a final remark, notice that making the choice I = t l / z in Eq. (jlj), the 
growth law takes the standard scaling form 




2.2 Two-time Quantities 

Once the role of A(e) is understood in the growth law, the next step is to 
establish how it affects the other observables in the problem. 

In this paper, we shall be concerned with spin systems which are quenched 
below the critical temperature at t = 0. For what follows, the basic two-time 
quantities are the autocorrelation function and the autoresponse function. 
The first of these is defined by 

c(t,t w ) = (^(O^)) - (^00X^(0). (7) 

Here, o~i = ±1 is the spin variable at the lattice site i, (t,t w ) are a pair of 
times after the quench, with the convention t w < t, and angular brackets 
(•) denote average over all sources of randomness in the problem: initial 
condition, thermal fluctuations and quenched disorder. After taking these 
averages C(t, t w ) is space-translationally invariant and there is no dependence 
on the choice of %. 

The autoresponse function is defined by 



5hi(t v 



(8) 

h=0 



where hi(t) is a space and time- dependent external field conjugate to <7j, 
and (-)h denotes the average in the presence of the field. This instantaneous 
response function is very difficult to measure, both numerically and experi- 
mentally. For this reason, time integration is used as a means to improve the 
signal to noise ratio [30] ■ We shall focus on the zero- field- cooled susceptibility 
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(ZFCS), obtained by switching on the external field from the time t w > up 
to the observation time t: 

X (t,t w ) = [ dt'R{t,t'). (9) 



The motivations and the advantages for using this form of integrated response 
function, with respect to others, have been discussed at length in Ref. [28] . 

Let us denote by 0(t, t w ) either one of the above observables. In the pure 
case and for quench temperatures low enough to neglect thermal fluctuations, 
the aging-scaling relation is obeyed jU EJ ED] 

0(t,t w ) = L~°<(t w )F (Jjj^Y (10) 

where a is the scaling exponent and F(x) is the scaling function. 

The extension of Eq. (fTOj) to the disordered case is a non-trivial problem. 
According to the SU hypothesis, disorder should affect only the growth law, 
leaving the form of Eq. (TXOT) unaltered with the same scaling function, that 
is 

0(t, t w , e) = L~ a (t w , e)F (j^y) • (H) 

In other words, once time is re-parametrized by L(t,e), there should be no 
difference between disordered and pure systems. This is a strong statement. 

Alternatively, one can allow for an explicit disorder-dependence of the 
scaling function 

0(t,t w ,e) = L-«p(^, -M, (12) 

and extend the crossover scenario to two-time observables. Here L(t, e) has 
been denoted by L, L(t w ,e) by L w and A(e) by A. Using the notation x = 
L/L w and y = X/L w , the x y region is controlled by the unstable y* = oo 
fixed point, while the x ^> y region is controlled by the stable y* = fixed 
point. Adopting the convention that hat-free symbols refer to the pure fixed 
point and hatted symbols to the disordered fixed point, if disorder is an 
irrelevant perturbation from the analysis presented in the previous subsection 
follows that the stable fixed point y* = is pure and that the general scaling 
relation in Eq. f|T2|) yields the limiting behaviors 

r-ocv \ I L~ a F(x), for x<^y, / 1QN 
L ™ F(x ' y ^{Ll°F( X ),ior X >>y, (13) 
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where F(x) ^ F(x). The exponent a is different from or equal to a, depend- 
ing on whether lim^oo F(x, y) is singular or not. Conversely, if disorder is a 
relevant perturbation the stable fixed point y* = is disorder-controlled and 
the limiting behaviors are reversed 

w 1 lV) ~ \ L~ a F(x), for x » y 1 } 

with a being different from or equal to a, depending on whether lim^o F(x, y) 
is singular or not. For instance, in the case of the autocorrelation function 
there are no singularities in F(x, y), since «c = &c — as a consequence of 
the compact nature of domains, both for pure and disordered systems. 

In the general framework outlined above, the SU hypothesis holds as an 
asymptotic statement when disorder is irrelevant, while SU is violated and 
the asymptotics of disordered and pure systems are different when disorder 
is relevant. 



3 Domain Growth in d = 1 Disordered Ising 
Systems 

Let us now discuss d — 1 disordered systems. These will provide clear realiza- 
tions of the possible scenarios outlined in the previous section. As a general 
comment, it should be kept in mind that even if d = 1 systems are param- 
agnetic at any finite temperature, as long as the domain size L(t) is smaller 
than the equilibrium correlation length £(T), the coarsening dynamics is the 
same as in the quench to T = where ferromagnetic ordering occurs. This 
is explained in detail in Refs. (101 HU [32] . 

We consider the ferromagnetic Ising chain with Hamiltonian 

N-l N 
i=l i=l 

The system evolves with non-conserved Glauber spin-flip dynamics [3]. In- 
terfaces in this case are point defects, separating domains of up spins from 
domains of down spins. In order to identify the length A, and to keep the dis- 
cussion simple, let us focus on the energetics of the particular configuration 



8 



[cr]j, containing a single defect breaking the bond Jf. 

Y°! (16) 
[ -1, for i > j. K ' 

The energy change when the defect moves from the bond j to another location 
r lattice spacings away with, say, r > 0, is given by 

AE jd+r = 2[J j+r -J j -Y / h i }. (17) 

i=j 

In the pure system, the ferromagnetic coupling is the same for all pairs of 
nearest neighbors (Jj = J > 0). In the absence of an external field (hi = 0), 
AEjj +r = for any r. The potential landscape is flat and defects perform 
free diffusion, yielding the behavior of Eq. ([1]) in the large-time regime [32J. 

In the d = 1 ferromagnetic RBIM, the couplings Jj are non-negative 
independent random variables, and hi = as in the pure case. For a given 
configuration of the couplings, there are now local minima and maxima in 
the potential landscape. A local minimum occurs when a defect breaks a 
bond flanked, to the right and to the left, by stronger bonds. Similarly, a 
defect sits on a local maximum when the opposite bond configuration occurs. 
The energy difference in this case is given by 

AEjj +r = 2(Jj +r — Jj), (18) 

which depends on the actual values of the couplings involved, but not on 
their distance. For instance, taking the Jj to be uniformly distributed in the 
interval (1 — e, 1 + e) with < e < 1, the upper bound on AEj d+r is given 
by 4e. Using the Arrhenius-type relation between escape time and barrier 
height £g ~ e EB ' T , there remains defined a characteristic time 

r ~ e 4e/T , (19) 

as the time needed to overcome the highest energy barrier, and a character- 
istic length scale 



A 



4e/(zT) _ 1 



(20) 



We take the Boltzmann constant ks — 1 throughout the paper. 

In the d — 1 RFIM, the ferromagnetic coupling is the same for all pairs 
of nearest neighbors (Jj = J > 0), as in the pure case, while hi = ±e is an 
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uncorrelated random field with expectation^! (hi) = and (hihj) = e 2 Sij. In 
this case, domain walls perform random walks in a random potential of the 
Sinai type [33]. The landscape contains minima and maxima, as in the d — 1 
RBIM. However, the crucial difference is that the height of energy barriers 
is not bounded and scales with the distance traveled by the defect: 

AE jjj+r ^--£hi-er^. (21) 

i=j 

The characteristic disorder length is obtained by balancing thermal energy 
with the barrier height 

A = (T/e) 2 . (22) 

Therefore, from Eqs. (120]) and (|22|) . it follows that when e — > the two 
possible scenarios, originating from Eq. ([5]) and analyzed in detail in Sec. |2j 
find realization in the d = 1 RBIM and RFIM, respectively. Indeed, this is 
what has been found [TUl [TT] . In the d = 1 RBIM disorder is irrelevant since 
from Eq. (120]) follows lim^o A = 0. Our numerical simulations [11] have 
shown very precisely that the growth law, the autocorrelation function and 
the ZFCS scale according to Eqs. §6§, (fT2l and (fTBl . with a crossover from 
pre-asymptotic disordered to asymptotic pure-like behavior. Actually, in the 
pre-asymptotic regime, the growth law is algebraic with a disorder- dependent 
growth exponent, indicating that A should more properly be regarded as a 
marginal scaling field. A similar behavior will be encountered in the d = 2 
RBIM discussed in Sec. HI Furthermore, the ZFCS exponent vanishes: 

a x = (23) 

for all disorder values, as in the pure case [32]. Hence, the ZFCS scaling 
function, like that of the autocorrelation function, does not have singularities 
for y — > oo. 

On the other hand, the d = 1 RFIM [10] is a realization of the opposite 
case, since Eq. (122]) shows that lim^o A = oo and, therefore, that disorder 
is relevant. The growth law, the autocorrelation function and ZFCS scale 
according to Eqs. (Q, f fT2]) and f ]T4]) and the crossover occurs from pure-like 
to disordered behavior. Again, as in the d = 1 RBIM, the exponents of 
the autocorrelation function and ZFCS vanish in the pure and disordered 

2 Clearly, here the average is only over the quenched randomness. 
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systems, revealing the absence of singularities in the scaling functions as 

Because of its importance in this paper, it is useful to recall the result for 
the growth law in the d — 1 RFIM. The limiting forms of the scaling function 
in Eq. are given by [10] 

C(y) ~ { ^7 )2 ' f ° r V « (24) 
yyj \ D, for y > 1. v 1 

Thus, in the pre-asymptotic regime (where y ^> 1 or L(t) <C A), the thermal 
energy exceeds the barrier height and growth takes place via free diffusion, as 
in the pure system (there is a possible dependence of the diffusion constant 
D on disorder.) In the asymptotic regime (where ?/ C 1 or L(t) ^> A), 
diffusion is of the Sinai type with Lit) ~ (lnt) 2 [33]. This is a consequence 
of the Arrhenius law with a barrier height which depends on L according to 
Eq. (12 ip. Then, the overall growth law, which interpolates between limiting 
behaviors and accounts for Eq. flUj) , can be written as 

f ~ L * e (W a . (25) 



4 Domain Growth in the d = 2 Random Bond 
Ising Model 

Let us next consider the d = 2 RBIM in the absence of an external field, 
with Hamiltonian 

W = -^JijViVj, (26) 

(ij) 

where the sum runs over nearest-neighbor pairs. The ferromagnetic couplings 
Jij are independent random variables with the same statistics as in the d — 1 
case, i.e., uniformly distributed in the interval (1 — e, 1 + e) with < e < 1. 

Let us first summarize the preexisting theoretical and numerical situation. 
According to the phenomenological HH [TS] theory, energy barriers in this 
system (when d > 1) ought to scale as a power of the domain size: 

E B (L) = kL^ } (27) 

where k is a disorder-dependent prefactor such that 

limK = (28) 
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and ip is an exponent dependent on dimensionality with the value ip = 1/4 
for d = 2. 

Let us, then, see what are the consequences on the basis of the gen- 
eral considerations made above. Identifying the characteristic disorder scale 
through the matching of thermal energy with barrier height, one has 

A = (T/k) 1 ^ (29) 

which implies lim^o A = oo. Therefore, in the HH theory disorder is relevant, 
which implies that the scenario found in the d — 1 RFIM should be replicated 
with a diffusive pre-asymptotic regime followed by asymptotic logarithmic 
growth 

(Int)^, for L»A (30) 
and with an interpolating formula analogous to Eq. fl25|) 

t ~ L z e (L/X)*_ (31) 



However, this HH scenario has never been actually observed. Rather, 
from a number of experiments on random magnets [32] and from extensive 
simulations of the d = 2 RBIM [T^l [2H [22] , the observed growth seems to be 
well-described by an algebraic law 

L(t) ~ (32) 

with a disorder- dependent exponent z. This growth law applies across the 
accessible time region, and there is no hint of a crossover. 

A possible explanation of the numerical findings has been put forward 
by Paul et al. [TP] on the basis of a logarithmic dependence of Eg on L, in 
place of the power law in Eq. fT27|) . However, Cugliandolo et al. [231121] have 
argued that the observation of a behavior consistent with Eq. (132]) does not 
necessarily rule out a growth law obeying Eqs. fl30|) and fl3T|) . since a disorder- 
dependent algebraic growth is compatible with Eq. ( 13"T|) as an intermediate 
regime, as it had been previously remarked by Bouchaud et al. [25] in the 
context of spin glasses. Thus, if we rewrite Eq. ([3T]) as 

I>~I* e {i<M+, (33) 

we obtain 

1 L* 

(34) 
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and, if L^/lnL can be treated as a constant c over the time interval of 
interest, the effective exponent is given by 

C CK /n , 

z - z+ x^ = z + T r' ^ ^ 

The present understanding of the observed algebraic growth ( 132]) is to a sort 
of standstill between these two mutually excluding explanations. 

In order to settle the issue, we have performed the most comprehensive 
simulations, to date, of coarsening in the d = 2 Glauber-RBIM. Our simula- 
tions have been done on N 2 square lattices up to maximum time t m Monte 
Carlo steps (MCS). Taking a random configuration of up and down spins, 
which mimics the disordered T = oo state before the quench, as initial con- 
dition for a run, we have let the system evolve with the updating rule which 
aligns spins with probability 1 with the Weiss field, if there is a majority of 
nearest neighbours pointing in the same direction, while the usual Glauber 
rule is used when the Weiss field is zero. This is equivalent to take couplings 
of the form = J+5ij, with 5ij uniformly distributed in (— e, e) and then let- 
ting T — > and e — > 0, while keeping e/T finite. Specifically, we have studied 
the coarsening dynamics with values e/T = 0,0.5,1,1.5,2,2.5. Kinetics in 
this form is appropriate for deep quenches, where it provides a considerable 
gain in the speed of computation, since updating trials are restricted to spins 
not aligned with the Weiss field, whose number roughly scales like 1/L(t). 
In the above limit disorder and temperature enter the transition rates in the 
form e/T and, therefore, also results depend on disorder and temperature 
only through the ratio e/T. Finally, we have checked that there is no dif- 
ference with the growth law and the aging properties of two-time quantities 
computed with the standard Monte Carlo algorithm. 

All statistical quantities have been obtained as an average over N run in- 
dependent runs. For each run, the system has a different initial condition 
and disorder configuration. The values of N, t m and N Tun for different values 
of e/T are listed in Table [TJ Notice that we have taken huge system sizes for 
e/T = and e/T = 0.5. This is necessary to obtain a reasonable time-window 
of domain growth before we encounter finite-size effects, as coarsening is more 
rapid for pure and near-pure systems with nonconserved dynamics. We have 
performed several checks to ensure that the data presented here are free of 
finite-size effects. This is crucial as the crossover phenomena we are inves- 
tigating are subtle and are expected to occur at very late times. We also 
stress that the t m - values for the present data sets are an order of magnitude 
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Disorder values 


System size 


Maximum time 


Number of runs 


e/T 


N 2 


t m (MCS) 


1 > run 


0.0 


12000 2 


1.01027 x 10 5 


20 


0.5 


8192 2 


3.12973 x 10 5 


20 


1.0 


4096 2 


1.006735 x 10 6 


75 


1.5 


4096 2 


1.006735 x 10 6 


75 


2.0 


4096 2 


1.006735 x 10 6 


75 


2.5 


4096 2 


1.006735 x 10 s 


75 



Table 1: Numerical parameters for d = 2 RBIM simulations for various 
disorder strengths. 

(or more) larger than those of Paul et al. 1 19 . 

Our findings for the growth law, the autocorrelation function and the 
ZFCS are presented in the following subsections. 

4.1 Growth Law 

We have obtained L(t) from the inverse density of defects, which is measured 
by dividing the number of sites with at least one oppositely-aligned neighbor 
by the total number of sitea^]. The plot of L vs. t in Fig. [1] shows the existence 
of two time regimes, separated by a microscopic time to of order 1. In the 
early time regime, for t < t , there is no dependence on disorder and growth 
is fast. The defects seeded by the initial condition execute rapid ballistic 
motion toward the nearby local minima. For t > to, the data sets in Fig. [1] 
appear to confirm previous observations [191 I2D [22] of a late-time regime 
with algebraic growth as in Eq. (132|) . Measuring slopes, the values of the 
disorder-dependent exponent z are displayed in Table [2] and Fig. [2} with an 
e/T dependence of the type observed in Ref. [22\. 

The improvement with respect to previous work comes from the quality 
of our data, which allow for a more refined analysis. In addition to estimating 
z we have extracted the effective exponent z eS (t) defined by 

1 - rf(lnL) (36) 



z eS (t) d(\nt) 



3 We have checked that this definition is in good agreement with the usual method of 
measuring L(t) by the half-peak-width of the equal-time structure factor. 
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io 3 



o—o e/T=0 




Figure 1: Growth law in the pure case (top curve) and for various disorder 
strengths. 



e/T 


z 





2 


0.5 


2.46 


1.0 


3.38 


1.5 


4.50 


2.0 


5.81 


2.5 


7.35 



Table 2: Exponent z for various disorder strengths. 
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Figure 2: Disorder-dependent exponent z vs. e/T. 
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In Fig. [5J we have plotted z e g(t) vs. L(t) and here comes the main point in 
the paper, since this plot allows to discriminate between algebraic and loga- 
rithmic asymptotic growth in favour of the latter choice. Indeed, if growth 
was described correctly by Eq. fl32l) . then z e s(L) ought to be independent of 
L, for L > Lq = L(to), with z e g(L) = z. It is clear from Fig. [3] that this is 
the case only in the pure system (e/T = 0). In all other cases, there is a slow 
time-increase which is not consistent with Eq. 032]) • In order to make this 
feature clear, in Fig. [5J we have drawn the straight horizontal broken lines 
corresponding to the values of z from Table [21 From the data follows i) that 
we can exclude the asymptotic validity of the algebraic growth law (|32|) (and 
with it of the logarithmic dependence on L of energy barriers) and ii) that 
the rather smooth and near-linear behavior of z e g(L), for the smaller values 
of e/T and L > Lq, is of the form 

z eS = C + hi", (37) 

where (,b,ip are fit parameters. Hence, inserting Eq. ( )3T|) into Eq. ( 156]) and 
integrating over time from i onward, we derive the growth law in the scaling 
form 

t h(L/X) 



where 

and where we have set 



t h(L /xy (38) 

h(x) = x c exp(x v ) (39) 
A = (<p/b)^. (40) 



We stress that in the above derivation there are no assumptions, other than 
the form ( 1571) of z c ff- Therefore, from the analysis of the numerical data we 
predict the existence of a crossover from algebraic to logarithmic growth, 
with the crossover length given by Eq. (T40]) . 

With the assumption Lo < A (to be verified a posteriori), and using 
h(x) ~ for x < 1, Eq. (158]) becomes 

{L/Xfe^ v = (L /X)<t/t (41) 

with the limiting behaviours 

Lit) - I Dtl/C > f ° r L K< A ' (42) 
L[t) ~ \ Apii(t/to)] 1/v , for L » A, l42j 
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Figure 3: Effective exponent z eS vs. L. Straight dashed lines indicate z. 
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e/T 


A 


c 


0.5 


762.2 


2.36 


1.0 


154.1 


3.15 


1.5 


64.5 


4.11 


2.0 


36.4 


5.13 



Table 3: Values of A and £ for various disorder strengths. 

where 

D = %■ («) 

Comparing with Eqs.( l30|) and ( 13T1) and identifying ip with ip, we find that 
the observed growth is of the HH-type, except for the replacement of the 
pure growth exponent z with the disorder- dependent quantity ( (which is 
distinct from z). The top line of Eq. (14*2]) accounts for the prior observations 
pill [2TJ [22] °f disorder-dependent power-law growth, without need of invoking 
the argument from Eq. (1331) up to Eq. (135]) . This is the pre-asymptotic regime 
followed by the asymptotically logarithmic one in the second line of Eq. (I42p . 
In addition, comparing Eq. ( 12 9 j) with Eq. ( I40p we can make the identification 

k — — — . (44) 

The first important conclusion, from the analysis of the numerical data 
presented above, is the validation, at least at the qualitative level, of the HH 
scenario. Clearly, it would be important to check also on the quantitative 
HH prediction ip = 1/4. However, it turns out to be difficult to obtain the 
numerical value of p from fitting the data sets in Fig. [3J We have considered 
three trial values p = 1/4, 1/2, 1, and for each one of them we have fitted 
(, b. In Fig. HI we have superposed these fits on the data sets for z e s vs. 
L for e/T = 0.5, 1, 1.5. In the following we use (p — 1, since this provides 
the best fit using the least-squares criterion. However, differences are too 
small to commit to one particular value and we leave the issue open. A 
precise determination of <p, as it is evident from Fig. HI requires to push the 
accuracy of the numerical computation far beyond what we have done. The 
corresponding values of ( and A = 1/6 are listed in Tableland plotted in 
the left and right panels of Fig. [51 Notice that since Lq ~ 1 , as it is evident 
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Figure 4: Comparison of z e g computed from Eq. fl37j) with ip = 1/4,1/2,1 
(continuous lines) and data from the simulations. 
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Figure 5: Exponent ( (left panel) and A (right panel) vs. e/T for if) = 1. 



from Fig. [H the assumed inequality L <C A is satisfied. Furthermore, A is 
found to decrease with increasing disorder as in Eq. (j22J) with 

A ~ (T/e) 2 . (45) 

Hence, comparing with Eq. f[2"§j) . we obtain k ~ e 2 /T. 

Summarizing, our data for z e g are evidence for the existence of a crossover 
of which, however, we can access only the pre-asymptotic region, since the 
values of the crossover length A are at the edge of the larger values of L(t) 
reached in the simulations for all disorder values (see Fig. [6]). Hence, numer- 
ically we can see only the onset of the crossover. Nonetheless, a number of 
conclusions can be drawn. As Eq. (|4"ST) shows, disorder is relevant implying 
the crossover to disorder-dominated behavior. The pre-asymptotic algebraic 
growth with a disorder-dependent exponent £ signals that disorder, although 
being globally relevant, acts marginally in the neighborhood of the pure (un- 
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T 



□ e/T=0.5 
« 6/T=l 




Figure 6: Fully developed crossover obtained extrapolating the growth law 
from Eq. (|38|) (continuous lines) with parameters ( and A from Table [3] and 
■0 = 1. Symbols depict the numerical data from the simulation. Asterisks 
mark the values of the crossover length A. 
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stable) fixed point Jj As remarked previously, this happens also in the d = 1 
RBIM, where a disorder-dependent exponent arises in the pre-asymptotic 
regime, but the physical mechanism producing these preasymptotic disorder- 
dependent exponents remains unclear. Finally, inserting the fitted values of 
the parameters into Eq. (138]) . the behavior of the growth law can be ana- 
lytically extrapolated for arbitrary large time. This has been done in Fig. [H] 
where, in addition to finding excellent agreement in the region where the sim- 
ulation data are available, the overview of the predicted crossover behavior 
is displayed. 

4.2 Autocorrelation Function 

Having established the existence of a crossover in the growth law, and that 
disorder is a relevant perturbation, we can anticipate that SU will not hold 
(see Sec. [2]). If the autocorrelation function did obey SU as in Eq. ffHJ|) . 
keeping in mind that ac = for pure and disordered systems, the plot 
of C(t,t w ) against L(t)/L(t w ) (for different values of e/T and t w ) should 
collapse all the data sets on a single master curve. Figure [TJ shows that this 
is not the case. Instead, there is collapse only when t w is varied keeping e/T 
fixed, but the master curves of these partial collapses, corresponding to the 
different values of e/T, are distinct. What we have called partial collapse has 
also been reported in the context of the diluted Ising ferromagnet jTTJ l3Tj . 
without reaching, however, a definite conclusion on the SU validity. 

The behavior in Fig. [7J can be accounted for on the basis of the two- 
variable scaling C(t,t w ,e) = F c (x,y) in Eq. f Ti"2"|) with ac = 0. The pure 
case F c (x, oo) = F c (x) corresponds to the lower black curve, where the 
collapse is very good since all the dependence on t w takes place through x. 
However, in the presence of disorder, y is finite. The remaining curves in the 
plot are explained considering that with the t^-values used in the simulations, 
the variation of y is small while keeping e/T fixed. This produces the quasi- 
collapse on the coloured curves, each one corresponding to a distinct value 
of e/T. But the variation of y becomes sizable when e/T (and thereby A) is 
changed, producing the separation of the curves with different colors. The 
trend in Fig. [7J shows that F c (x,y) increases as y decreases, for fixed x. 
Hence, the bundle of curves generated by letting y vary is bounded below 

4 The statements that the perturbation is marginal and that the pure fixed point is 
unstable are compatible when nonlinear terms in the expansion around the fixed point are 
taken into account [36 . 
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Figure 7: Autocorrelation function for different disorder values. The t w values 
have been taken according to the rule t w (k) = t m i n (t m / 2t m , n ) ( fc ~ 1 )/ 10 where 
^min — 200, k is an integer running from 1 to 11 and t m is the maximum 
time in Table [TJ Accordingly, the waiting times are logarithmically spaced 
between t w = 200 and t w = t m /2. 
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by the pure scaling function Fc(x), and must be bounded above by the 
numerically inaccessible scaling function Fc(x) of the disorder-dominated 
fixed point. Exactly the same behavior is observed in the autocorrelation 
function of the d = 1 RFIM [10]. 

Within the framework of our findings, we conclude that previously re- 
ported observations of SU in the autocorrelation function for the d = 3 
RFIM [H [T5] and the d = 2 RBIM [22] are due to time-ranges and disor- 
der values deep inside the pre-asymptotic regime, where it is not possible to 
detect the onset of the crossover. 

4.3 Autoresponse Function 

The qualitative difference between the autoresponse and autocorrelation func- 
tions is that the exponent a x in the scaling relation 

X {t,t V) ,e) = L- a *F x (x,y), (46) 

is not constrained to vanish, like ac, by geometrical requirements. Hence, 
in principle, the crossover could take the most general form of Eq. (|T4|) with 

Before analysing the data it is useful to summarize the structure of scaling 
in the pure case [2B], where Eq. (146"j) reads 

X (t,t w ) = L- a *F x (x) (47) 

and for large x the scaling function decays with the same exponent appearing 
in the prefactor 

F x (x)~x~ ax . (48) 

As a matter of fact, the measurement of the decaying slope of x( x ->tw) for 
large x is the fast numerical method to extract the scaling exponent a x . The 
physical origin of this feature can be understood introducing the effective 
response associated with a single interfac^], which is defined by [27] 

X(t,t w ) =p(t)xes (t,t w ) (49) 

where p(t) = L^ 1 ^) is the interface density. Then, from the scaling form in 
Eq. (|47D follows 

XeS (t,t w ) ~ L£- a ^xF x (x). (50) 
5 For instance, a configuration with a single interface has been introduced in Eq. (|T6|) . 
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On the other hand, if there is only one interface in the system the response is 
not expected to age, yielding a time-translation invariant Xcs(t — t w ). There- 
fore, regarding the left hand side as dependent on time trough L(t — t w ) and 
considering the large x regime where it is possible to replace L(t — t w ) by 
L(t), the dependence on L w must drop out of the right hand side implying 
xF x (x) ~ x( 1-Q! x). This, in turn, implies the validity of Eq. ( I48p and 

Xe s(t,t w ) -L 1 ^*. (51) 

Let us now go to the data in presence of disorder. We have computed 
x(t, t w , e) with the algorithm derived in Ref. [38], which does not require the 
switching on of the external perturbation, thus improving the precision and 
the efficiency of the numerical computation. The use of algorithms of this 
type is indispensable for the study of linear response when the system is 
quenched to T = 0. The quantity L^ x x(t,t w ,e) has been plotted in Fig. [8] 
against L/L w for different values of e/T and t w , using for all disorder values 
a x = 0.6, which is the value producing the best data collapse in the pure 
case (top black curves). The structure of the plot follows the same pattern 
as in Fig. [7J with good partial collapse when e/T is kept fixed and distinct 
master curves when e/T is varied, for the same reasons explained in the 
discussion of the autocorrelation function. The collapse of the curves with 
the pure exponent a x shows clearly that the numerical data belong to the 
preasymptotic region controlled by the pure fixed point up to the onset of the 
crossover, since the y dependence of the scaling function F(x,y) is manifest 
in Fig. E 

Furthermore, since the different curves in the plot appear to be obtained 
one from the other by a vertical shift, in this region of parameters the de- 
pendence on x and y can be taken to be of the multiplicative type 

F(x,y) = f(x)g(y) (52) 

with f(x) decaying algebraically for large x, exactly as in Eq. (jlSl) for the pure 
case. In order to highlight this feature, in the inset of Fig. [5] we have plotted 
F(x,y)/F(x = 2,y) obtaining the superposition of the different curves, that 
is the elimination of the y-dependence, which is a nice check on the validity 
of the assumption made in Eq. fl52|) . Therefore, the curve in the inset is 
the normalized plot of f{x), whose dacay in agreement with Eq. ( I48p is a 
confirmation i) that scaling is controlled by the pure fixed point and ii) of 
the validity of the single interface structure of the response. 
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Figure 8: Zero-field-cooled susceptibility for different disorder values and t, 
values chosen according to the rule in the caption of Fig. [7J 
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Figure 9: Zero-field-cooled susceptibility for a single interface vs. L(t), taken 
from the data of Fig. [H for different disorder values. Asterisks mark the 
values of the crossover length A. 



As a further check of this latter point, we have measured the ZFCS 
Xsingie {t,t w , e) after preparing an initial state with a single interfac^] and 
for different disorder values. The plot in Fig. [9j which also includes the pure 
case, does indeed show that Xswg\e(t,t w ,e) grows with a power law in agree- 
ment with Eq. (|5Tl . The exponent does not seem to depend on e/T and it is 
consistent with the expected value 1 — a x ~ 0.4. Conversely, in front of the 
power there is a prefactor A dependent on e/T (see Table HJ), which enters 
into the factor g(y), in Eq. (152]) . and is responsible for the separation of the 
master curves in Fig. [Sj 

6 The initial single interface is prepared along a diagonal in the lattice, in such a way 
that spins on the interface experience zero Weiss field. If the interface was prepared along 
the horizontal or vertical axis, the Weiss field would not be zero and no flip could take 
place with the accelerated algorithm presented at the beginning of this section. 
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e/T 


A 


0.0 


0.00235 


0.5 


0.00184 


1.0 


0.00125 


1.5 


0.00093 



Table 4: Amplitude A of Xsingie(^> t w ) extracted from the data of Fig. for 
different disorder values. 

In summary, since our data do not allow to access the full crossover 
region, we cannot make statements on the interesting question whether the 
exponent a x will be different or not from a x when the disorder-dominated 
asymptotics are reached. However, precisely because what we observe is 
the pre-asymptotic regime controlled by the pure fixed point, this result is 
important in itself, as explained in detail in Ref . [TT] , since it brings new and 
independent evidence on the value of a x in the pure system. In fact, together 
with the d — 1 result in Eq. (123]) . this result substantiates the conjecture put 
forward [27J for the pure system that a x depends on dimensionality according 
to 

a x = (d-l)/2, (53) 

when d < 3. The numerical finding a x ~ 0.6 does not exactly coincide with 
the value a x = 0.5, predicted by Eq. fl53|) for d — 2. This means that either 
the phenomenological formula or the precision of the numerical computa- 
tion still need improvement. However, it seems clear that the alternative 
conjecture [23 001 ED] predicting 

a x = 1, (54) 

with no dependence on dimensionality [ID], is ruled out by the data from 
the d — 1 and d = 2 RBIM. Preliminary results on this point have been 
presented in Ref. [TT] . 

5 Summary and Discussion 

Let us conclude this paper with a summary and discussion of the results 
presented here. We have reported results from the most comprehensive sim- 
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ulations to date of domain growth in the d = 2 ferromagnetic RBIM with non- 
conserved Glauber kinetics. We undertook this computationally expensive 
study with a specific purpose in mind, viz., to obtain a clear understanding 
of possible crossovers in the domain growth law and other statistical prop- 
erties of the evolution morphology. As regards the domain growth law, our 
results are incompatible with algebraic growth scenario of Ref. [H] and are 
compatible with a crossover from a pre-asymptotic regime having power-law 
growth with a disorder-dependent exponent to an asymptotic logarithmic 
regime [L(t) ~ (lni) 1 /^]. The pre-asymptotic regime is consistent with ex- 
isting results for the RBIM [191 [22] • To the best of our knowledge, the 
predictive power of our result for the late-stage result is new and allow us 
to give numerical support to the HH scenatio. However, in spite of the huge 
numerical effort, we cannot conclusively determine the exponent (p which 
characterizes the logarithmic regime. Clearly, more work is needed in this 
direction. 

We use RG arguments to understand the crossover in the growth law 
in terms of the competition between a pure fixed point and a disordered 
fixed point. Within this framework the corresponding crossover in two-time 
quantities, like the autocorrelation function C(t,t w ,e) and the autoresponse 
function t w , e), is expected. The scaling function for C(t, t w , e) shows the 
onset of the crossover from a pure form to a disordered one. Even though 
the numerical data are limited to the pre-asymptotic region, there is clear 
evidence that the SU hypothesis does not hold for two-time quantities. A 
similar behavior is observed in the data for x(t,t„,,e), where the crossover 
region is not entered deep enough for F x (x,y) to develop the y-singularity 
which eventually should change the scaling exponent from the pure value a x 
to the disorder value a x . Nonetheless, the measurement of the pure exponent 
a x ~ 0.6 is quite an interesting result in itself, adding new evidence in support 
of the phenomenological conjecture of Eq. (1531) . 

At a more general level, the RG framework formulated in this paper, in 
the context of the d = 2 RBIM, proves to be comprehensive enough to give 
also a natural explanation for the known results for the d = 1 RFIM and 
d = 1 RBIM. As such, it is a good candidate for the understanding of coars- 
ening in disordered systems, whose predictive power could be tested at higher 
dimensionality or in other systems, like the site diluted Ising model, where 
the crossover region perhaps is more easily accessible than in the RBIM [17j. 
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